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, Abstract. We show, using a covariant and gauge - invariant charged multifiuid 

perturbation scheme, that velocity perturbations of the matter - dominated dust 
Friedmann - Lemaitre - Robertson - Walker (FLRW) model can lead to the generation 
of cosmic magnetic fields. Moreover, using cosmic microwave background (CMB) 
constraints, it is argued that these fields can reach strengths of between 10 -28 and 
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CO , 10 29 G at the time the dynamo mechanism sets in, making them plausible seed field 
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candidates. 
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1. Introduction 

Cosmological magnetic fields have generated considerable debate ever since they were 
first observed [1]. Such large scale fields occur in galaxy clusters, high redshift 
condensations, spiral and disc galaxies, and have strengths between 1CU 7 to 1CU 5 G [2]. 
While the structure of the magnetic fields in spiral galaxies appear to indicate that they 
were generated and sustained by a dynamo mechanism [3], the origin of the necessary 
seed field is still unknown. For the dynamo mechanism to be successful, the present day 
strength of these seed fields must lie within a range of 10 -23 [4] to 10~ 30 G [5], the lower 
value being for a dark energy dominated spatially flat universe. 

Primordial magnetic fields are of particular interest since they explain both 
the fields seen in nearby galaxies as well as those detected in galaxy clusters and 
highly redshifted condensations. Several different schemes generating these fields 
have been suggested, many of them based on breaking the conformal invariance of 
electromagnetism. Since this can be achieved in more than one way, this explains the 
variety of the proposed mechanisms in the literature, such as coupling the photon to 
a scalar field either during inflation or in the subsequent era of preheating [6] (see [7] 
for a critique), or assuming the breakdown of Lorentz invariance either in the context 
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of string theory and non - commutative varying speed of light theories, or due to the 
dynamics of large extra dimensions [8]. Although these proposals successfully explain 
the generation of primordial fields, it is often done at the expense of simplicity because 
they extend the physical parameter space. 

Recently it has been suggested that if a magnetic seed field exists below the 
necessary strength for the dynamo mechanism to work, gravitational waves could be 
used to pump the magnetic field by many orders of magnitude [9]. The mechanism uses 
the primordial gravitational wave spectrum generated during inflation, and couples this 
via Maxwell's equations to the magnetic seed field. This model has the attractive feature 
that one need not introduce new physics or extend the number of fields considered. 
Of course, a disadvantage is that it relies on the yet unobserved inflationary, model - 
dependent gravitational wave spectrum, and that still requires another model to create 
the weak seed field. 

With this in mind, we propose a mechanism based solely on the physics of self - 
gravitating plasmas. Since velocity and density perturbations naturally occur in the 
early Universe, it is interesting to examine whether such perturbations can induce 
magnetic fields which can be sustained at appreciable levels until the onset of non - linear 
gravitational collapse (for a discussion of density perturbations in the presence of weak 
magnetic fields see, e.g., [10-12]. The mechanism is similar to Harrison's protogalaxy 
model [13], and the Biermann battery effect [14], in the sense of yielding vorticity driven 
magnetic fields, but we note that the battery effect in our formalism would be of second 
order, while the Harrison effect relies on Thomson scattering. 

Using a recently developed formalism for describing perturbations in a general - 
relativistic multi - component plasma [15], it is found that the value of the magnetic 
field generated by primordial velocity perturbations is well within the limits required if 
the dynamo mechanism is to work. The limitations of the method, as well as the results, 
are then discussed. 

2. Preliminaries 

To begin with, we present some of the key equations from [15] which were 
obtained by linearizing the exact Einstein - Maxwell equations using a two - parameter 
approximation scheme characterized by 

• e g - gravitational: a ab , u ab , Vo- ab ,Vu ab , etc. 

• e em - electromagnetic: E a ,B a . 

In this way, terms which are second order in the gravitational variables, e 2 ,, second order 
in the induced electromagnetic fields, e 2 em and all cross terms e g e em are neglected. 

2.1. First Order Multifluid equations 

Using the above approximation scheme, the Einstein - Maxwell equations are linearized 
around a pressure - free FLRW model. Hence all spatial gradients and velocity 
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components perpendicular to the fundamental observer with 4-velocity u a must vanish 
in the background. This implies that inhomogeneities and velocity perturbations play a 
first order role. In addition, the isotropy and homogeneity of the background demands 
that the electromagnetic field and accordingly the total charge density p vanishes in the 
background. 

From Maxwell's equations, we find the first order equations 

E {a) - curl5 a = -e l j {a) - §6£ a , (1) 

B {a) + curl£ a = -§OB a , (2) 

V a E a = e - V , (3) 

V a B a = , (4) 

where curl_B a = e abc VbB c and 9 is the cosmological expansion defined through the scale 
factor S: = 3S/S 

Neglecting all but electromagnetic interactions between the fluids (e.g. collisions), 
then, for each fluid species (i), conservation of energy, momentum and particle density 
leads to the following equations: 

£(*) + A*(*)V U(i) = -fi(i)Q , (5) 

ft® + n (i) V ' a v a {{) = -©n (i) . (7) 

Since we are working in the energy frame, where the total heat flux vanishes, the equation 
for total energy conservation is simply given by 

p = -Qp. (8) 

Because no thermal effects are included, we may express the fluids' energy density as 
Pi = min,i and the energy frame choice further implies y"^t>° = 0. 

i 

Specializing to two components, it is convenient to introduce the following variables 
N = m + n 2 , (9) 
n =m-n 2 , (10) 

v = |K+t£) , (ii) 

v a =\{vt-v a 2 ) . (12) 
The following system of first order equations then follows, using equations (6) and (7): 
n = -NV a v a - &n , (13) 

N = N—V a v a - QN , (14) 
p 

u<°> = -- (— + — 1 E a - -Qv a , (15) 

where we defined Ap — p\ — p 2 and the total energy density p — p\ + p 2 , respectively. 
Since we have pi = m^i, the latter can be written as 

p — |(mi + m 2 )N + \{m 2 — mi)n, (16) 
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a relation which will be used very often in the following analysis. The Raychaudhuri 
equation then becomes 

e = -|0 2 -i(m 1 + m 2 )N - \{ mi - m 2 )n . (17) 

It is useful to introduce the following quantities, namely 

y - £, (is) 

« 2 = TT^—, (19) 

6eomim2 

P 2 = —r^ v (20) 

£0(^1+^2) 

The fraction Y = n/N obeys a propagation equation, 

y + fey + \c?ilY = 0, (21) 

which may be derived from the above definitions and equations (13)— (17). 
3. Velocity Induced Electromagnetic Fields 

For a cold plasma, the currents for each fluid species may be written as 

where is the charge and vfy is the velocity of the species under consideration. Since 
we require the plasma to be neutral on the whole, the species are of opposite charge. 
Hence, the total current j a appearing in Maxwell's equations reads to first order 

3 a = fi + ti = -eNv a . (23) 

From Maxwell's equations (1) -(4), using (23) and (15), one can then deduce second 
order wave equations for the induced electromagnetic fields. They are 

E {a) - V 2 E a + l&E {a) + [§6 2 + (fa 2 + |) E a = 2(3 2 fi (v a Y - |0t; a ) , (24) 

B {a) - V 2 B a + |eS {a) + [|6 2 + i/i] B a = -2(3 2 fi curl v a . (25) 

Observe that B a and curlf a are both purely solenoidal, whereas V a Y has no solenoidal 
part. It is worthwhile to note that the magnetic field is solely sourced by inhomogeneities 
in the velocity in contrast to the electric field which is sourced by inhomogeneities in 
the number density and velocity perturbations. Both equations look strikingly similar, 
the differences originating either from the total current or from a gradient in the charge 
density (in the case of V a F). The additional 3/4a 2 -term in the electric wave equation 
comes from the non-stationarity of the total current and its largeness - a 2 ~ 10 42 for 
an e + e~ - plasma - leads directly to the high-frequency behaviour of plasma effects, as 
will be shown below (see also [15] for a discussion of the high-frequency plasma mode 
in the gravitational instability picture). 

It will be useful to introduce expansion normalized variables, 

&a = -q, SSa = -Q-, JCa = Q ■ (26) 
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r + 3 )/i 



e 



£ a = 2[3 2 £ [ V a Y - ^Qv a ) (27) 



1„ 3a 2 A . 

v {a) + -ev a = -__e<f a , 



Ha) 



Ha) 



ie 2 -i/i)^ a = -2/3 2 /iX, 



3a 2 



3 20 



(28) 
(29) 
(30) 



Equation (30) follows from (26) using (28) and Maxwell's equation (2). In order to 
find solutions to the above equations, we extract from them the scalar and solenoidal 
(vector) parts (cf. Appendix B) and solve them separately. 

3.1. Scalar modes 

In analogy with (B.3), we set V = SV a v a and S = SV a £ a . Equation (28) then 
transforms into 

v + \ev = -~e* = |aV?r; (3i) 

where the last equality is a direct consequence of Maxwell's equation (3). Combining 
Y = —V/S with (31) and using (8) and (17) together with 

SV a V 2 <& = VV+ (-§0 2 + (32) 

one can show that the scalar part of the electric wave equation (27) reduces to 



e+ I |e-| ]<? + 



<f = 0. 



(33) 



It is also easy to see that equation (31) additionally gives rise to propagation equations 
for V and Y: 



9 V4 



6^ 



v = o, 



Y + -GY + - A a 2 ^Y = 0. 
3 4 



(34) 
(35) 



Hence, equations (33)-(35) all emanate from (31). We now specialise our considerations 
to a flat FLRW background with zero cosmological constant, for which \x = 1/3 2 and 
= 2 ft always hold, solutions may easily be obtained: 



V(r) 
«f(r) 
Y(t) 



= -j= jV; cos(cj lnr) + ^ + V?J sin(u; lnr)| , 



--tL—i(2Vi + 3V ) cos a; lnr + ^ l - sin u; lnr 

9 CT y'T I Old 

U 1 1 f, rtT , otwx / , x ( 2 - 18a 2 )^ + 3V? . , , , 
^T^TTTe i^K + 3 l/J cos(a; lnr) + 

q u sin (a; lnr J 



(36) 
(37) 
(38) 



Cosmic magnetic fields from velocity perturbations in the early Universe 



6 



Here, we introduced the dimensionless time-coordinate r = t/U, where tj denotes some 
arbitrary initial time. The scale factor becomes now S(r) = S^r 2 / 3 . Initial conditions 
of the velocity perturbation are chosen to be Vi = V(l) and V- = V'(l) (a prime stands 
for d T ). The frequency of the solutions is proportional to uo = a/« 2 — 1/36 and grows 
logarithmically in time. The solutions show the same high - frequency behaviour that 
was obtained in [15]. 



3.2. Vector modes 



According to (B.6), we set S a 

"s ka) - v 2 i a + (e - 1) i {a 



Scm\S a etc., and obtain from equations (27)-(30) 



+ 



1 3 a 2 



1(a) 



- v 2 ^ a + 



(a) 



3 26 



a " 8^ 



-e 2 + 
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3 2 1 



3 29 



/i 



(39) 
(40) 
(41) 
(42) 



Specialising to a flat FLRW background and performing a harmonic decomposition (see 
Appendix B), these equations become 



S" + — <f ' + 
3r 



r 4/3 

3a 2 



3r 
4 

+ 37 



4(3 2 t 
L 2 



+ 



S 



a 



+ 



JtT' + — JT 
3r 



,-4/3' 

3a^ 

8/3 2 



f + 



3r 2 
1 



^ 2 



9r 2 



3r 



(43) 

(44) 

(45) 
(46) 



where we have dropped the index V [denoting that the variables in equations (43)-(46) 
are vector harmonic coefficients], and we have defined r = tjU and used 

2 /a \ \2 -t -, 



kt{ 
~S 



47T 

T 



Xh 
A 



.4/3 



.4/3 



(47) 



for the contribution of the Laplacian in (39) and (41), respectively. In (47), \ H = l/H 
is the Hubble length, A is the physical wavelength associated with the comoving 
wavenumber k = 2nS/\ and the index i stands for evaluation at initial time ij. If 
the wavelength of the mode is much greater than the horizon, eg. Aj ^> Aj^, we may 
neglect the terms containing L 2 . However, we can neglect that term in (43) and the 
system (45) - (46) throughout, because the a 2 - term dominates as long as rL 3 <C a 3 
holds and this criterion is only violated for very late times or ultra - short wavelengths. 
It follows that the above equations can then be solved analytically and the general (real) 
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solutions are found to be 

v (t) =Ci + ^{C 2 cos(cjlnr) + C 3 sin(cjlnr)}, (48) 

^( r ) = -l^i - {(^ + 6cuC 3 ) cos(culnr) + (C 3 - 6coC 2 ) sin(wlnr)} , (49) 

9 or 9 cr 

jT(r) = -^{Dicos^lnr) + Z^sin^lnr)} , (50) 
1 8 Z? 2 1 

<^(r) =/J 3 - T7 ^ + -^- T7 ^{/J 1 cos(^lnr)+ J D 2 sin(^lnr)}, (51) 

where u = \J a 2 — 1/36 and = a/ a 2 /2 — 1/36. Since we think of the electromagnetic 
fields as being generated by the velocity perturbations, we choose initial conditions for 
v and JtT, respectively: namely, Vi = v(l), v[ = v'(l) and J% = J^' = J^'(l). 

Notice that the system of equations (43) and (44) is equivalent to a third order ODE for 
v (or $), for any value of L, while the system of equations (45) and (46) is equivalent 
to a second order ODE for , for L = only. However, the vector modes of Jf a are 
linearly related to those of v a [according to (26)], thus we have SQJt ~ v. Therefore 
consistency requires the vanishing of the integration constant C\. For this set of initial 
conditions, the integration constants Cj and D { become 

C 2 = v „ C 3 =Hi±M, D ^, Zl , D ^'B±ML, (52) 
Zuj boo 

while there is no restriction for D 3 , since we put L ~ 0, but with the assumption that 

the initial magnetic field vanishes, we get a relation between D% and J%. The solutions 

for the velocity perturbation and the (expansion normalized) electric field agree then 

with those found in the scalar case and show the same behaviour in time, as expected. 

The induced expansion normalized magnetic field attains two parts, a standard decaying 

part and a weakly decaying oscillatory part due to the plasma. 

4. Applications 

4-1- Generated fields 

Observe that the magnetic field (51) is rather slowly decaying and therefore still could 
play a role in some astrophysical processes under favourable conditions. If we assume 
that the velocity induced magnetic field vanishes initially, then we may approximate its 
magnitude by the following expression: 

W^^^-^e)^ (53) 

where is the magnitude of the initial velocity curl perturbation. Restoring SI units, 
we find for the physical magnetic field 

1/6" 



S X / N / L \ 1/6 

"Ml \ I m l \ I m e \ ^ I H 



\B(t)\< ^ )( — ) e )jtr il 

- 1 mj \m e J \m 1 + m 2 J \t 



ix2xlO" 7 G, (54) 
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Figure 1. Generated magnetic field in units of 10 27 G, for the initial and final rcdshift 
1000 and 100, respectively. 



or alternatively, using 1 + z — (1 + Zj)/r 2 / 3 , 

\B(z)\<jeihi±^) i - ' 1 



l + Zi 



1 + Zi 



;i + zf' 2 x 1(T 24 G, (55) 



where we also have employed = 3H = SH (1 + zf' 2 with H = /i(9.8Gyrs) _1 , and 
neglected the mass factor. f Thus, velocity curl perturbations of magnitude J% ~ 10~ 5 
in an e + e~ - plasma starting after decoupling, when Thomson scattering becomes 
negligible [z ~ 1000), would induce a magnetic field with a strength between 10~ 26 
and 10~ 27 G at a redshift of z ~ 100 (see figure 1 above), a redshift well within the 
limits before nonlinear effects become important. Redshifting to z ~ 10, at the onset of 
the dynamo mechanism [1] reduces the field strength to between 10~ 28 and 10 -29 G.§ 

The above argument is also applicable for the decaying electric field (for both scalar 
and vector modes). If we require the velocity induced electric field to vanish initially, 
then its envelope is given by 

Af3 2 1 

KnvO - )! = -— , (56) 

o a 

where i>j is the magnitude of the initial scalar velocity perturbation. We resort once 
again to SI units and find for the physical field 

|£„vWI =(rSL)(nH)( -=M * (j) " 2 7 x 2 x l<r 3 Vm", (57) 
\m e J \m e J \m 1 +m 2 J \t J t 

or alternatively, 

\E env (z)\ = Vih { ^ Z J /4 x 2 x 10- 23 Vm- 1 . (58) 

| The mass parameter in (55) and (57) is always of the order of one for an electron-positron- or an 
electron-ion-plasma. 

§ The constraint concerning the initial velocity perturbation and Vi, respectively, stems from the 
standard CMB results [16]. 
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Hence, using the same data as above, i.e., velocity perturbations of magnitude Vi ~ 10~ 5 
in an e + e~-plasma starting after decoupling [z ~ 1000), would lead to an induced 
electric field of strength ~ 10~ 26 Vm _1 at z = 100. 

Comparing the energy density of the induced electric and magnetic fields, we 
find that E 2 nv /c 2 B 2 ~ 10~ 15 , thus showing the well known fact that the electric field 
contribution, due to Debye shielding, is negligible in a cosmological context. 

5. Summary and discussion 

In this paper, we have considered the induced electromagnetic fields due to velocity 
perturbations in a charged multifluid, using a gauge - invariant, covariant approach. 
We have, in a self - consistent manner, investigated the behaviour of the interacting 
fluids and electromagnetic fields in the case of a flat FLRW background. In the matter 
dominated era, due to the non - vanishing vorticity of the velocity perturbations in the 
two - component fluid, we find that there is a net magnetic field with a magnitude 
between 10~ 30 and 10 -31 G today. Since velocity perturbations are naturally occurring 
in the early Universe, this field thus represents a suitable candidate for a seed field, which 
could be amplified by the galactic dynamo. Moreover, the model is self - consistent, and 
does not invoke any other physics than general relativity and classical electrodynamics 
The notion of cosmic magnetic fields is nowadays generally accepted, and 
observations seems to indicate their presence on all scales of the Universe. On the 
other hand, the existing standard models for the generation of these fields [1] all require 
a seed field, and it is still somewhat of a mystery as to where this primordial field stems 
from. There have through the years been a number of suggestions for the origin of the 
seeds, many of them making use of yet to be confirmed physics. At the same time, CMB 
data is getting more and more detailed [17], and we now have a very good handle on 
the size of the different types of perturbation that occur in the early Universe. Since 
it is highly plausible that the plasma state is a good approximation of the cosmological 
fluid at certain stages of the evolution of the Universe, the model presented here lends 
itself naturally to the analysis of a possible source of the much sought-after magnetic 
seed field. In many of the classical mechanisms for generating the galactic and extra- 
galactic magnetic fields, the vorticity of the fluid plays a crucial role in the generation 
of the magnetic fields [13, 14]. This, in conjunction with the aforementioned CMB 
data, gives our model the attractive feature of being not only internally consistent, but 
also consistent with present day cosmological data. On the other hand, scattering and 
thermal effects have been neglected, and would, if included, surely contribute to the 
dynamics of the fluids in a nontrivial and important way, but these issues are left for 
future studies. We note however that including a radiation gas and Thomson scattering 
would bring the model to a form more accurately describing the plasma at pertinent 
epochs, as well as bringing it closer to Harrison's protogalaxy model. 
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Appendix A. Some commutator expressions 



In this section, we give some useful expressions for commuting spatial derivatives, up to 
first order, in the case of dust spacetimes (see also [18]). We will assume in the following 
that ViV and X a are first order quantities. 

V 2 iv) ' = V 2 A - ^0V 2 A , (A.l) 

(A.2) 

(A.3) 

(A.4) 



v [( ,v,, v,a- = ( l -(-y 2 ^// j • . 



V a V 2 A = V 2 V a A + -6 2 - -fi) V a A , 



. 9 

v [a v fe] v c x d = ( -ie 2 + \im ) h c [a v b] x d + h d [a v c x b 



v a v 2 x b = v 2 v a x b + -e 2 - -n 
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V a X b + V b X a - h ab V c X c 



Se abc V b V 2 X c = V 2 (Se abc V b X c ) 



abc\ 



(A.5) 
(A.6) 



Appendix B. Isolating scalar and vector modes 

We define our harmonics as eigenfunctions of the Laplace-Beltrami operator [19,20] 

v 2 g = -^g, q = o, (b.i) 

where Q stands for a scalar, vector or tensor harmonic. For example, a first order vector 
field X a may be expanded covariantly in scalar and vector harmonics 

X a = X s Q a s + X v Q a v , (B.2) 

where an implicit summation in this expansion is understood. In order to extract the 
purely scalar modes, X s , of a first order vector field, one basically takes the divergence 
(multiplied by the scale factor S for convenience) and readily obtains the following 
relations: 

X = SV b X b = X s (kQs), (B.3) 
X = SV b X b = X s (kQ s ), (B.4) 
X = SV b X b = X s {kQ s ). (B.5) 

The solenoidal modes, X v , can be obtained by applying curl and noting that the 
curly harmonics Qy = Se abc V b Q^ also satisfy relation (B.I). The relations analogous 
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to the scalar case are now 




(B.6) 
(B.7) 
(B.8) 



Thus, for fixed comoving wave number k, X and X s as well as X a and X v will obey 
identical equations. We like to stress that all relations above are valid within the limits 
of our two-parameter approximation scheme. 
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